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ABSTRACT 


A  semi -infinite  plate  is  exposed  to  the  stepwise-decelerated 
hypersonic  flow  at  zero  angle  of  attack.  The  boundary  layer  reconfigu¬ 
ration  gives  rise  to  an  induced  pressure.  To  obtain  this  temporal  pressure 
distribution,  the  transformed  momentum  and  energy  equations  have  been 
solved  using  the  approach  of  Ref.  [3]  where  the  solutions  are  obtained  in 
the  range  0<t<1  .  These  solutions,  with  extension  to  the  range  1  <t<]  . 3 , 
are  used  to  evaluate  the  displacement  thickness  for  the  unsteady  boundary 
layer.  The  differentiated  displacement  thickness  is  finally  used  in  the 
induced  pressure  expression  involving  the  acoustic  approximation.  A  family 
of  curves  is  obtained,  for  the  initial  zero  pressure  qradient  case,  cover¬ 
ing  the  ranqe  0<g  <1  +  C^-)  M  ^  ,  5<M  <30  and  0<x<1.3.  The  shearinq 
stress  and  heat  transfer  at  the  wall  are  obtained  for  the  indicated  range 

of  q  and  M  and  0<t<1.8.  The  absolute  maqnitude  of  the  transient  contri- 
"W  e  - 

bution  to  the  induced  pressure  distribution  is  found  to  be  oreater  for  the 
case  of  no  heat  transfer  at  the  wall  and  0<q  <1,  and  smaller  for  q  >1  as 
compared  to  the  corresponding  contribution  for  the  case  of  stepwise-acceler¬ 
ated  flow  treated  in  Ref.  [3].  Some  attempts  have  also  been  made  to  improve 
the  solution  of  the  momentum  equation  by  modifying  the  numerical  method  of 
Ref.  [3]. 
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CHAPTER  I 


STATEMENT  OF  THE  PROBLEM  AND  ASSOCIATED  LITERATURE 

1 . 1  The  Unsteady  Hypersonic  Flow  Problem 

The  theory  of  unsteady  compressible  boundary  layers  has  not  been 
extensively  studied  in  the  past.  One  of  the  reasons  for  this  is  that, 
in  many  engineering  problems,  boundary  layer  growth  occurs  in  so  short 
a  time  that  it  is  often  sufficient  to  assume  the  steady  state  although 
the  speed  changes  with  space  and  time.  There  are  cases,  however,  where 
a  more  detailed  investigation  is  desirable.  These  cases  include,  for 
instance,  oscillating  wings,  blades  rotating  in  non-uniform  air  streams, 
unsteady  nozzle  flow  and  aircrafts  and  missiles  in  unsteady  flight.  The 
maneuverabil ity  and  control  of  aircrafts  and  missiles  is  completely  de¬ 
pendent  upon  the  forces  acting  on  the  components  of  these  aerodynamic 
vehicles.  In  evaluating  these  forces  viscous  interactions  must  be  con¬ 
sidered.  The  research  on  the  blunt-leading-edge  and  the  sharp-leading- 
edge  bodies  is  the  best  source  of  information  on  viscous  interaction 
phenomenon . 

The  problem  to  be  considered  herein  is  the  interaction  between 
the  viscous  and  inviscid  effects.  The  initial  zero  angle  of  attack,  high 
Mach  Number  flow  is  stepwise-decelerated  by  a  small  amount  in  the  direction 
of  the  initial  flow.  In  order  to  obtain  the  interaction  induced  pressure, 
the  development  of  the  velocity  and  enthalpy  boundary  layers  will  be  con¬ 
sidered  first.  From  this  growth  of  the  boundary  layers,  equations  will  be 


1 


2 


derived  which  will  account  for  the  effects  of  the  interaction  between  the 
viscous  and  inviscid  flow  upon  surface  skin  friction,  heat  transfer  and 
the  normal  forces  on  the  plate. 

Within  a  hypersonic  boundary  layer,  the  temperature  is  expected 
to  be  high,  which  means  that  the  density  will  be  low  and  the  equation 

6*  -  /“(i  -  T^n-)  dy  (1-1) 

0  pe  e 


serves  to  indicate  that  the  displacement  thickness  will  be  large.  Ordi¬ 
narily,  this  would  be  of  no  particular  consequence  except  close  to  the 
leading  edge  of  a  sharp-nosed  object.  Near  the  sharp  leading  edge  where 
the  displacement  thickness  is  increasing  from  zero,  the  turning  of  the 
flow  external  to  the  boundary  layer  causes  an  initial  compression  accompa¬ 
nied  by  a  shock  wave.  This  compression  with  the  following  flow  expansion, 
in  turn  affects  the  build-up  of  the  boundary  layer,  and  thus  we  have  a 
problem  of  interaction  between  the  external  inviscid  flow  and  the  boundary 
layer. 


In  the  hot  gas  layer  near  the  surface  it  is  certainly  not  at  all 

clear  at  first  whether  the  usual  Prandtl  Boundary  Layer  Equations  are  appli- 

2 

cable  for  hypersonic  flow.  However,  Shen  [1]  has  shown  that  when  (6/x)  «1 , 

the  higher  order  viscous  and  heat  conduction  terms  [(yu  )  ,  (kT  )  ]  in  the 

x-momentum  and  energy  equations  are  small  compared  to  the  terms  (yu  )  and 

y  y 

(kT  )  .  In  addition  the  static  pressure  ratio  across  the  boundary  layer 

»y  y 

is  negligible. 


The  problem,  then,  is  to  determine  the  solutions  for  a  classical 


•  i 


Prandtl  boundary  layer  subjected  to  an  "external  pressure  gradient"  which 
is  not  known  in  advance,  but  which  depends  on  the  growth  of  the  boundary 
layer  itself. 

Lees  and  Probstein  [2]  have  divided  the  total  flow  region  into 
two  main  regions:  I,  a  "leading  edge"  region;  II,  a  region  of  "no-slip" 
as  shown  in  Fig.  1.1 


EDGE 

Fig.  1.1 

Hypersonic  Flow  Over  a  Flat  Plate 

Breaking  the  leading  edge  region  into  two  sub-regions  indicated 
by  1  and  2  in  Fig.  1.1,  they  point  out  a  relatively  small  domain  "1"  around 
the  leading  edge  which  is  of  the  order  of  20  mean  free  paths  in  breadth 
and  domain  "2",  the  slip  flow  region. 

In  region  II,  as  pointed  out  by  Lees  and  Probstein,  slip  and 
temperature  jump  at  the  surface  are  negligible,  and  the  continuum-flow- 
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concepts  apply  throughout.  The  no-slip  region  is  further  subdivided  into 
three  domains:  (3)  strong  interaction  region;  (4)  transition  region; 

(5)  weak  interaction  region.  Based  upon  the  observations  of  the  authors 
of  Ref.  [2],  the  strong  interaction  region  should  be  important  only  for 
flight  Mach  Numbers  greater  than  about  15. 

In  the  present  analysis  we  shall  confine  ourselves  to  the  weak 
interaction  regime.  As  pointed  out  previously  the  development  of  the 
unsteady  boundary  layer,  especially  in  hypersonic  flow,  has  received  little 
attention.  The  review  of  the  literature  related  to  the  problem  under  con¬ 
sideration  is  presented  in  the  following  section.  Next  the  two-dimensional 
nonsteady  comoressible  boundary  layer  equations  are  transformed,  using 
the  Dorodni tsyn-Howarth  transformation  and  the  stream  function.  The  use 
of  a  similarity  transformation  reduces  the  number  of  variables  from  three 
to  two.  The  resulting  momentum  and  the  energy  equations  are  coupled  and 
can  not  be  solved  by  ordinary  numerical  techniques.  These  equations  may 
be  linearized  keeping  in  view  the  small  perturbation  nature  of  the  problem. 
The  linearization  technique  used  here  is  similar  to  that  of  Ref.  [3], 

For  zero  pressure  gradient  case  the  linearized  momentum  and  the  energy 
equations  become  uncoupled  but  still  the  nature  of  the  momentum  equation 
does  not  permit  a  closed  form  solution.  Therefore  the  numerical  method 
for  its  solution  was  adopted.  Ordinary  type  of  numerical  methods  did 
not  succeed.  A  summary  of  these  methods  has  been  included  in  the  appendices. 
The  iterative  finite  difference  method  of  Appendix  'C1  was  found  to  give 
the  desired  solution  of  the  momentum  equation.  However,  the  method  was 
too  costly  computer  time-wise.  Instead  the  numerical  method  developed  in 
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Ref.  [3]  was  used  with  some  changes  in  the  computational  technique.  With 
these  changes  the  solution  of  the  momentum  equation  was  available  for 
0<t<1 . 3 .  The  solution  of  the  energy  equation  was  obtained  for  the  adia¬ 
batic  wall  as  well  as  for  finite  heat  transfer  at  the  wall  from  the  existing 
literature  [3]. 


The  shearing  stress  at  the  wall  was  directly  obtained  from  the 


solution  of  the  momentum  equation.  The  differentiated  displacement  thick¬ 
ness  was  used  to  obtain  the  normal  pressure  distribution  involving  the 
acoustic  approximation.  The  variation  of  the  steady  state  boundary  layer 
thickness  was  also  obtained  for  zero  and  non-zero  heat  transfer  at  the 
wall.  The  relevant  results  have  been  discussed  in  the  last  chapter. 

1 .2  Review  of  Previous  Work 

Blasius  was  the  first  (see  paragraph  65,  Ref.  [21])  to  investigate 
the  unsteady  laminar  boundary  layer  flow  over  various  bodies  starting  from 
rest.  Several  other  investigators  also  analysed  the  unsteady  boundary 
layer  flow  for  special  cases.  F.K.  Moore  [22]  included  in  his  treatment 
the  unsteady  boundary  layer  flow  for  longer  times  elapsed  since  starting 
and  for  higher  speeds.  He  considered  the  case  of  compressible  laminar  flow 
over  a  semi -finite  flat  plate  in  rectilinear  accelerated  flight  through 
still  air.  In  his  analysis  he  found  that  a  group  of  parameters 
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can  be  formed  incorporating  the  coordinate  along  the  surface  and  the  free 
stream  velocity  and  its  derivatives.  For  uniform  acceleration  (Ug  =  AT), 
the  quantities  (1.2)  reduce  to  a  single  paramater  c  =  x/AT  (=-).  For 
C»1 ,  the  effect  of  the  leading  edge  has  not  been  felt  at  the  station  x 
and  the  classical  "starting-from-rest"  solution  applies.  If  c«l ,  the 
growth  of  the  boundary  layer  with  x  must  be  considered  and  Ref.  [22] 
points  out  that  the  flow  may  be  regarded  as  nearly  quasi-steady;  that  is, 
at  any  instant  the  motion  is  nearly  that  which  would  be  obtained  in  steady 
flow  at  the  conditions  prevailing  at  that  instant.  Moore  considered  an 
insulated  flat  plate  problem.  Ostrach's  [6]  extension  includes  the  heat 
transfer  over  an  isothermal  surface. 

K.  Stewartson  [7]  studied  the  fluid  motion  induced  by  the  uni¬ 
form  motion  of  a  semi-infinite  flat  plate  in  its  own  plane.  Using  the 
stream  function  defined  by 

i|>  =  Mvt)^2  <|>(c,t)  (1.3) 


with 


Stewartson  expands  <j>  in  a  power  series  in  t  with  coefficients  that 
are  functions  of  c  .  This  series  has  as  the  leading  term  the  error 
function  solution  of  the  Rayleigh  Problem  in  the  sense  of  ||-  =  erf  j  , 
t-*0,  and  all  subsequent  terms  are  identically  zero.  He  then  indicates 
that  an  essential  singularity  is  likely  to  exist  at  t=1  ,  £-*».  However, 


.  ...  J 
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no  such  singularity  may  be  expected  for  0<t<1.  He  also  points  out  that 
the  Rayleigh  solution  valid  for  t<1  ,  may  be  joined  nonanalytical ly  at 
t=1  to  some  x-dependent  solution  valid  for  t>1 .  Stewartson  further 
explains  why  the  x-dependence  can  not  enter  via  a  power  series  to  describe 
the  transition  from  a  function  of  (y,t)  in  Uoot<x  to  a  function  of 
(y,t,x)  in  U0ot>x.  In  order  to  explain  the  disappearance  of  t-dependence 
as  t-*»  where  <f>  must  approach  the  Blasius  function,  he  mentions  that  the 


terms  depending  on  t  may  be  exponentially  small. 


Fig.  1.2 

Stewartson' s  Impulsive  Plate  Problem 


I 

•'  • 


. 
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In  his  later  publication  [8]  Stewartson  shows  that  the  asymptotic 
behavior  for  large  t  (Blasius)  is  not  relevant  to  the  solution  in  t<1 
(Rayleigh).  Considering as  the  dependent  variable  and  taking  the  coef- 
ficient  of  the  mixed  derivative  as 


he  finds,  with  the  help  of  (1.3),  that  the  momentum  equation  in  the  x- 
direction  has  the  character  of  a  heat  conduction  equation  in  t<1  and  <J> 
is  independent  of  x  and  given  by 


||  =  erf  (§)  ;  t<1  .  (1.4) 

Physically  this  means  that  the  fluid  at  any  point  (x,y)  of  the  boundary 
layer  is  not  "aware"  of  the  existence  of  the  leading  edge  until  t  =  x/U^  , 
the  time  taken  for  a  signal  to  travel  from  the  leading  edge  with  the  main¬ 
stream  velocity  relative  to  the  plate.  But  when  t>1  ,  then  k  ,  regarded 
as  a  function  of  c  is  partly  positive  and  partly  negative.  Hence  small 
disturbances  now  travel  in  both  directions  (increasing  and  decreasing  x) 
so  that  the  solution  must  depend  on  conditions  at  x=°°  and  at  x=  1 . 

Again  we  may  say  that  the  fluid  at  (x,y)  when  t>x/Uoo  is  "aware"  of  the 
leading  edge  and  its  previous  (downstream)  history. 

The  work  of  Lam  and  Crocco  [9]  gives  details  of  the  growth  of 
boundary  layer  induced  by  a  shock  travelling  down  a  semi -infinite  flat 


* 
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plate.  At  t=0  the  shock  is  at  the  leading  edge  of  the  plate,  x=0. 

For  the  special  case  of  a  very  weak  shock  the  problem  becomes  identical 
to  Stewartson's  impulse  plate  problem.  Using  the  Prandtl  boundary  layer 
equations  under  the  Crocco  transformation  they  have  given  sufficient 
conditions  for  unique  solution  to  the  governing  equations  for  the  case 
py  =  const.  The  treatment  of  the  equations  is  consistent  with  the  physics 
of  the  problem  corresponding  to  situations  in  a  shock  tube.  The  boundary 
layer  is  divided  into  two  parts:  I,  a  region  where  the  old  solution  of 
Mirels  [10]  is  valid  (t<1);  II,  the  region,  t>1 ,  where  the  leading  edge 
effects  and  the  effect  of  the  shock  solution  downstream  both  must  be  con¬ 
sidered.  t  is,  again,  here,  taken  as  Ugt/x  . 

Cheng  [11],  while  considering  the  problem  of  a  sharp-edged  flat 
plate  started  from  rest  in  an  incompressible  fluid,  expands  the  stream 
function  in  powers  of  the  quantity. 


he  obtains  an  approximate  solution  downstream  of  the  leading  edge  and 
points  out  that  the  solution  can  be  continued  as  far  downstream  as  one 
wishes.  The  formal  difficulty  of  the  appearance  of  "an  essential  singu¬ 
larity"  at  £=1  and  n=°°  ,  encountered  by  him  and  Stewartson  in  the 
method  of  perturbing  the  Rayleigh  type  solution,  disappears.  Cheng  in  his 
later  publication  along  with  Elliot  [13]  included  in  his  treatment  an 
arbitrary  plate  velocity  and  the  effects  of  compressibility. 


10 


Schuh's  [23]  approach  is  based  on  the  use  of  approximate  mo¬ 
mentum  integral  method  (Pohl hausen-Method  in  the  wider  sense).  Using  the 
parameter 


T  = 


he  obtains  the  growth  of  the  boundary  layer  over  a  circular  cylinder  and 
also  the  development  of  the  boundary  layer  and  the  frictional  drag  over  a 
flat  plate  of  finite  length  for  the  case  of  impulsive  start  from  rest  and 
for  the  case  of  steady  acceleration.  In  the  expression  for  x,  Uq  is  the 
constant  reference  speed  and  c  is  the  characteristic  length.  Schuh  noticed 
that  there  is  a  discontinuity  in  slope  of  the  skin  friction  distribution 
for  a  flat  plate  in  the  range  l<x<°°. 

On  examination  of  the  literature  pertinent  to  the  unsteady  com¬ 
pressible  boundary  layers,  it  becomes  evident  that  x  should  be  one  of 
the  parameters  to  be  used  in  such  analysis.  It  is  indicated  in  reference 
[3]  that 


n  =  /  [ 


(m+1)  Ue 
2vref  cx 


]  Y. 


Uet 

and  x  =  — 


are  the  proper  variables  to  transform  the  non-steady  boundary  layer  equations 
to  similarity  forms. 

Based  on  the  "Piston  Theory"  of  Lighthill  [14],  which  was  sup¬ 
plemented  by  Landhal  [15]  and  Miles  [16],  it  has  been  pointed  out  in  Ref. 


11 


[3]  that 


au  ah 

e  _  _ e  . 

at  "at 


(1.5) 


Giving  details  of  the  above  expressions  it  has  been  mentioned  that  one 
obtains  an  equivalent  of  the  Ackeret  pressure  distribution  instantaneously 
on  the  plate  at  large  Mach  Numbers. 


CHAPTER  II 

THE  GOVERNING  EQUATIONS 

2.1  General  Form  of  the  Equations 

The  Prandtl  boundary  layer  equations  in  the  unsteady  form  are 
taken  as  the  equations  governing  the  flow  in  a  compressible,  unsteady, 
laminar,  two-dimensional  boundary  layer.  They  are  [8,  17]: 


(2.1) 


(2.2) 


|£+|_(PU)  +|_(pV)  =0 


(2.4) 


with  the  boundary  conditions 


when  y  =  0, 


u  =  v  =  0 


} 


when  y  =  ye 


h  =  hw(x,t) 
u  =  Ue(x,t) 


} 


} 


(2.5) 


7  1 
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The  pressure  in  the  boundary  layer,  being  independent  of  y  , 
may  be  given  by  its  value  of  y=yg  according  to  the  inviscid  solution: 


at 


+  u  — 

e  ax  '  ax 


(2.6) 


where  Ug  denotes  the  prescribed  potential  motion. 

2.2  Dorodnitsyn-Howarth  Transformation 

Using  the  Dorodnitsyn-Howarth  transformation  [8] 


Y  =  Jy£-dy  (2.7) 

0  p°° 

and  defining  the  stream  function  as 


the  boundary  layer  equations  become 


a2^  ,  a^  a2i^ 
ayat  ay  axay 


M  iJt  =  Ps.  u  +  cv  ^ 
3x  sy2  p  ue  3x  3y3 


(2.9) 


ah  ,  a^  ah  a^  ah 
at  ay  ax  "  ax  ay 


pe  ai \>  1 1  8Ue  .  C  a2h 
p  ay  ue  3x  -  P r  9Y2 


+  VooC  (^|)2  (2*1°) 

ay 
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where 


(2.11) 


(2.12) 


It  has  been  assumed  that 


C  =  C(x) 


VP 

H»Pco 


(2.13) 


and  that  Pr  and  Cp  are  constant.  The  product  yp  is  evaluated  at  some 
point  within  the  boundary  layer  so  that  it  is  independent  of  Y  . 


2.3  Similarity  Forms 

The  following  forms  of  the  similarity  variables  [3]  may  be  used 
in  order  to  reduce  the  independent  variables  from  three  to  two. 


(m+1)  Up 
n  =  7  [2v  cx 


■]  Y 


(2.14) 


2v  cx  U 

*  =  /  [  ;-;+1)  e]  f(n.x)  (2.15) 


where 


uet 


T 


X 


(2.16) 
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and 


(2.17) 


In  addition,  let 


h  =  g(n,T)  G(x,t) 


(2.18) 


On  substituting  (2.14)  and  (2.15)  into  the  differential  equations  (2.9) 
and  (2.10) ,  one  obtains 


f'  !T+  f  -  Uef,fl  7+  f'f'T  ar-  fcff,'Ue  +  fUef"  J 


-  ff  "t 


aUe  1  9Ue  h  aUe  2  Ue 

8x  2  TT  3x  ti  3x  -  x  T  U 

e 


(2.19) 


and 


Ue  3G  Ue  Ue  3G 

+Tm_eGaf.  +g|iuef' 


,  U  ,  U  u  .  u 

i^-Gfg1  -  Im  ^  fg'  +  x  ^  fGg'  -  xfGg'm  ^ 


2  x 


fe  2  ft  Ue  (m+1 )  Ue  Ggl,  ,  (m+1)  Ue  f,,2 

'  p  ue  T  01  x  2Pr  x  bg  2  x  T 


(2.20) 


m+1 
2c  v 


00 


where 
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f  = 


9f 

8t 


f'  = 


3f 

3n 


and 


$  = 


9'  = 


l£ 

3t 


1SL 

3n 


Let 


U  =  ex 
e 


m 


(2.21) 


6  = 


2m 

m+1 


(2.22) 


and 


G(x,t)  =  G(x)  =  he(x) 


(2.23) 


Using  the  energy  equation  for  the  free  stream  and  expression  (2.6),  one 
obtains 


_ e  _  3 d  _  .I  e 

pe  3x  3x  peUe  3x 


and  employing  the  power  law  distribution  of  expression  (2.21),  one  may 


wri  te 
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3r=-mi-  <2-24> 

Substituting  (2.22)  through  (2.24)  into  (2.19)  and  (2.20)  one 
obtains,  under  the  assumption  of  perfect  gas,  the  equations  given  by 
Rodkiewicz  [3],  namely, 

(2-e)f‘  +  2(l-3)x(ff 1  1  -  f'f)  -  f'"-ff"  -  3(g-f'2)  =  0  (2.25) 

Pr(2-3) [1  -  20i||Tf.]  $  .  g- •  .  Pr[f  -  2(l-3)xf]  g' 

-  Pr(y-l)  Me2  f ' ,2  =  0  (2.26) 


where 

u  2  2 

(y-1)  Me2 
e 

For  the  case  of  weak  interactions  (zero  pressure  gradient)  equations  (2.25) 
and  (2.26)  reduce  to 

2f '  +  2x(ff"  -  f'f')  -  f'"  -  ff"  =  0  (2.27) 

2Pr  (1-Tf)  §-g"  -  Pr(f-2xf )g 1  -  Pr(y-l)  M02  f"2  =  0  (2.28) 
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2.4  Linearization  of  the  Problem 

The  flow  pattern  immediately  following  the  start  of  an  impulsive 
motion  in  a  viscous  fluid  is  necessarily  the  same  as  would  occur  in  an 
inviscid  fluid,  except  for  the  presence  of  (infinitely  thin)  vortex  sheets 
along  all  solid  boundaries.  This  follows  from  the  fact  that,  in  the  case 
of  a  suddenly  "created"  velocity,  the  dominant  terms  in  the  Navier-Stokes 
equations  are 


f^—lgradp  (2.29) 

so  that  curl  u  =  0,  i.e.,  a  potential  flow  pattern  would  be  found  except 
near  solid  boundaries  where  the  no  slip  condition  causes  the  initial  vortex 
sheet. 

An  interesting  generalization  of  this  result  was  pointed  out  by 
Sears  [18],  who  considered  an  impulsive  change  superimposed  on  an  already 
established  flow  pattern.  He  showed  that  the  initial  motion  following  any 
impulsive  change  of  the  boundary  conditions  consists  of  the  superposition 
of  the  velocity  pattern  existing  just  before  the  change  and  the  inviscid 
flow  velocity  pattern  due  to  the  impulsive  boundary  values  (together  with 
the  corresponding  infinitely  thin  wall  vortex  sheet).  This  is  the  result 
of  equation  (2.29)  which  dominates  the  pattern  at  the  moment  of  change. 

The  terms  of  this  equations  are  all  linear. 

This  would  mean  that  the  velocity  pattern  following  the  step- 
deceleration  of  a  flat  plate  moving  in  Hypersonic  viscous  flow  may  be  re- 
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presented  as: 


The  dotted  line  indicates  the  final  steady  state. 

Let  us  assume  that 

f(n,x)  =  f 2 (n )  -  Af(n,T)  (2.30) 

where  f 2(n)  satisfies  the  final  steady  state  equations  existing  in  the 
literature  [13,14]. 

This  type  of  linearization  technique  was  used  by  Rodkiewicz  [3] 
for  a  step-accelerated  flat  plate  in  Hypersonic  flow. 

Substitution  of  equation  (2.30)  into  equation  (2.27)  results  in 

2(Af ) 1  +  2x[f“(Af)  -  f£(Af)']  -  (Af)'" 

-  f2(4f)"  -  f£'(Af)  -  0  (2.31) 


where  the  second  order  terms  in  Af  have  been  dropped  and  the  following 


. vsof  X  nf  etafq  Jfi.Vbe*  3  i-<  ;i-2 
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steady  state  momentum  equation  has  been  used 

fZ  '  +  f2  fz  =  0  (2.32) 

Also,  the  steady  state  energy  equation  is 

g£'  +  Pr  f2  gj  +  Pr(Y-l)  Mg2  f"2  =  0  (2.33) 

2.5  Associated  Boundary  and  Time  Conditions 

When  y  =  0,  Y  =  0,  n  =  0,  u  =  0,  f'(0,T)  =  0,  f2(0)  =  0 
and 

( Af ) 1 (0 ,t)  =  0  (2.34) 

when  y  =  0,  Y  =  0,  n  =  0,  ip  =  0,  f(0,x)  =  0,  f2(0)  =  0  and 

Af (0 ,t)  =  0  (2.35) 

when  y  =  ye,  Y  =  Y0,  n  =  ne»  u  =  Ue2,  f'(ne.T)  =  1,  f2(ne)  =  1  and 

(Af)'  (ne,x)  =  0  (2.36) 

when  t  =  0,  t  =  0,  u  =  u-|  -  (U^  -  Ue2), 
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U  ,  U  , 

f'(n>0,0)  =  Jr1  -  (IT1  -  1)  and 
1  ue2  ue2 

U  U  , 

( Af ) 1  (n>0,0)=  -  Tr1+  (rr1  -  D  +  fo(n) 

1  ue2  ue2  ^ 


This  may  be  integrated  to  give 

u  i  u  l  u  i 

Af(n>0,0)  =  -  /[rp— ]  f o ( n )  +  irr —  1  )n  +  f«(n) 

ue2  c  ue2  ue2  1 

when  y=0,  Y=0,  n=0,  h=h  and 

w 

g(o,x)  =  7 p- 
ne2 

and  for  the  case  of  an  insulated  plate 

g'(0,T)  =  0 

when  y  =  ye  ,  Y  =  YgS  n  =  ne,  h  =  hg2  and 

g(ye,T)  =  l 

when  t  =  0,  t  =  0,  h  =  h-|  and 

hl 

g(n.O)  =  t- — 
ne2 


(2.37a) 


(2.37b) 


(2.38) 


(2.39) 


(2.40) 


(2.41) 


■- 
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At  this  stage  it  would  be  useful  to  note  the  relation  between 
the  quantities  belonging  to  initial  and  final  steady  state  conditions. 
From  the  relation  (2.14),  one  may  write 


n 


1 


v  /u 

00 


el 


(2.42a) 


n2  ■  ^  Y  ^e2 


(2.42b) 


From  above,  we  obtain 


U  2  U  « 

n0  =  hi  /  n —  and  dn9  =  /[n^-] 
"  1  Uel  1  Uel 


(2.43) 


If  f-|  and  fg  satisfy  equation  (2.32)  with  the  boundary 
conditions: 


when  n-|  =  0,  f^  =  0,  f-j  =0 

(2.44) 

when  nl  =  ne’  =  1 

when  =  0,  f^  =  0,  f£  =  0 

(2.45) 

when  n2  =  ne’  f2  =  1 


The  two  solutions  must  be  the  same  and  accordingly  f-|(n-|)  and  *^2^2) 
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may  be  related  by  the  relation 


Uo9 

=  'f:2(v/[y^]n1  ) 


If 


u 


el 


0.01 


Then, 

f]  (n-j )  =  "f 2 ( 0 •  995n ) 


(2.46) 


(2.47) 


(2.48) 


This  may  be  represented  diagrammatical ly  as: 


Similarly  the  velocity  profile  may  be  represented  as: 


Blasius  Function  f 1 (n) 
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From  the  figure  one  can  easily  see  that 


fjCn,)  =  fj(n2)  =  f2(0-995rii) 


(2.49) 


CHAPTER  III 


SOLUTIONS  TO  THE  GOVERNING  EQUATIONS 

3.1  Numerical  Solution  to  the  Momentum  Equation 

The  solution  to  equation  (2.31)  with  the  time  dependent  boundary 
conditions  (2.34)  through  (2.37)  for  the  case  of  a  stepwise-decelerated 
flat  plate  in  viscous  hypersonic  flow  is  presented.  In  reference  [3], 
where  the  case  of  a  stepwise-accelerated  flat  plate  has  been  considered, 
a  numerical  solution  to  equation  (2.31)  with  the  appropriate  boundary 
conditions  is  discussed  in  detail.  The  difference  in  the  two  problems 
(of  stepwise-accelerated  and  stepwise  decelerated  flat  plate)  arises  from 
the  boundary  condition  (2.37).  The  author  of  reference  [3]  has  suggested 
the  recurrance  relation 


(3.1) 


where 


-a 


R 


m,n 


m,n  -  Tb  +  c  R  +  d  R  R  TT 
*  v  m,n  m,n  m,n-l  m,n  m,n-l  m,n-2' 


(3.2) 


and 


(c  C  T+d  R  0  C  ,  +  d  C  0) 
*  m,n  m,n-l  m,n  m,n-2  m,n-l  m,n  m,n-2/ 

Tb  +  c  R  i  +  d  R  \  R  T) 

v  m,n  m,n  m,n-l  m,n  m,n-l  mjn-2' 


(3.3) 


and  where 


R  -i  =  C  -I 

m ,  I  m ,  I 


=  Cm  9  =  0,  Rm  9  =  0.25 
m,2  m,2 


(3.4) 
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from  the  boundary  conditions  (2.34)  and  (2.35). 


Vn  >  bm,n  •  cm,n  '  dm,n  and  em,n  are  defined  in  APPendix 


•c . 


For  n  =  N  ,  where  the  "e"  condition  prevails,  the  following 
relation  has  been  suggested  in  view  of  the  boundary  condition  (2.36). 


Tn,N-l  4-R 


'm,N-2 


m,N 


m  ,N-2 


4-R 


m,N-2 


-  Rm,N-l 


(3.5) 


The  computations  are  based  on  the  following  shape  of  the 

molecule: 


V  i 

)n  +  1 

V 

/ 

}  \ 
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V 

m-1  m 

- ►  I 

Fig.  3.1 

7-Point  Molecule 

However,  the  numerical  method  did  not  produce  consistent  stream 
function  increment  (Af)  distribution  for  x>l .  Some  jumps  were  noted  in 
the  value  Af  .  A  search  was  made  by  the  present  writer  to  find  out  the 


3ll,  tnq  ion  bfb  bodJw  ><W  .^v»  on 
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"blow-up"  causes.  It  was  found  that  the  term 

/pr - Rm  N-l 

H  m,N-2  m,IN  1 

in  expression  (3.5)  goes  through  zero  as  t  gets  greater  than  one. 

An  iterative  finite  difference  method  was  tried  to  get  rid  of 
the  previously  mentioned  difficulty  for  t>1 .  The  method  did  succeed 
in  giving  values  upto  t  =  1.2.  However,  the  numerical  method  was  too 
costly  computer  time-wise  and,  therefore,  it  was  given  up  finally.  A 
summary  of  this  method  is  included  in  the  Appendix  ‘ C ' . 

It  was  decided  to  use  the  numerical  method  of  Ref.  [3]  with  some 
improvement  in  the  computational  technique.  In  order  to  reduce  the  "Round¬ 
off"  errors,  double  precision  was  used  in  the  computations.  With  this 
change,  Af*  distribution  was  obtainable  up  to  t  =  1.3  and  the  shear 
stress  at  the  wall  could  be  computed  up  to  t  =  1.8. 

In  place  of  using  the  previously  indicated  7-point  molecule, 
some  attempts  were  made  to  use  the  following  11-point  molecule: 


*See  Appendix  ' F '  for  the  curves  of  Chapter  III 
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This  molecule  is  the  result  of  applying  the  commonly  used  ex¬ 
pressions  for  central  derivatives.  The  third  order  derivative  in  in¬ 
direction  is  based  now  on  points  (m,n+2),  (m,n+l),  (m,n-l)  and  (m,n-2) 
instead  of  (m,n+l),  (m,n),  (m,n-l)  and  (m,n-2).  The  first  order  derivative 
in  the  x  -  direction  incorporates  points  (m+l,n)  and  (m-l,n)  in  place  of 
(m,n)  and  (m-1  ,n)  and  the  mixed  derivative  involves  points  (m-l,n+l), 
(m+l,n+l),  (m+l,n-l)  and  (m-l,n-l)  in  place  of  points  (m,n+l),  (m,n-l) 
(m-l,n-l)  and  (m-l,n+l).  The  11-point  molecule  would  permit  laraer  steps 
in  the  x  -  direction. 

Three  different  numerical  approaches  were  made  for  the  11 -point 
molecule  computations.  The  first  involved  a  recurrence  relation.  The 
second  used  the  Gauss's  Elimination  method.  The  last  utilized  an  iterative 
finite-difference  method.  All  the  three  methods  failed  to  give  a  con¬ 
sistent  Af  distribution.  A  summary  of  all  these  methods  is  included  in 
the  Appendices  'A', *  1 B 1  and  1 D 1 . 

3.2  Solution  to  the  Energy  Equation 

The  energy  equation  (2.28)  is  the  same  as  the  one  used  for  an 
accelerated  flow  problem  of  Ref.  [3].  Consequently,  the  following  solutions, 
for  Pr  =  1 ,  may  be  used. 

(1)  For  the  case  of  no  heat  transfer  at  the  wall 

g  =  i  +  Ac^ll  Me2  (i-f'2) 


(3.7) 


93 p f c  nt  (n,  T-i  )  bns  n0^^vb  *  T  *rt*  rt*  | 
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(2)  For  the  case  of  non-zero  wall  heat  transfer 


g  =  1  +  -^^Me2  (f'-f'2)  +  (gw-l)(l-f) 


provided  that  the  wall  temperature  is  kept  constant. 


3.3  Shearing  Stress  at  the  Wall 

The  shear  stress  at  the  wall  is  given  by 


Tw  “  9y'y=0 


or 


where 


w 


■  [y  (4)2  4k._ 


9n 


2  V9Y'  9yJy=0 


Using  (2.7),  (2.14)  and  (2.15)  one  obtains 


ii  3/2 

Vw  Ue  , 

Tw  T~  /~[2v  ex]  f  (0’T) 

r'oo  u  oo  "* 


xw  =  a,  f"(0,x) 


ii  3/2 
Vw  Ue 

“l  p„  cxj 


A  plot  of  the  shear  stress  as  a  function  of  t 


(3.8) 


(3.9) 

(3.10) 

(3.11) 

is  indicated 


in  Fig.  3.8. 
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3.4  Heat  Transfer  at  the  Wall 

The  heat  transfer  at  the  wall  may  be  obtained  from 


or 


qw  =  ('k  3y) 1 y=0 


qw 


k  h  ,  .u 

w ,e  n'  (o  T)f—  —1 
C  9  'U'  M3Y  3y;  1  y=0 


(3.12) 

(3.13) 


Differentiating  (3.8),  one  obtains 

g'  =  Me2f"(l-2f‘)  -  (gw-l)f ' '  (3.14) 

Using  (2.7),  (2.14),  (2.15)  and  (3.14)  one  finds  from  (3.13) 


k  h  P 

q  =  -  /[ 


2v6  cx^  2 


(y-1  ) 


’-g  +i] 
»w  J 


or 


where 


qw  =  a2f"(0,i) 


.  ,k.Wh.ePW  yr _ J! _ ir  _1  _  (y-1  )  M  2-, 

a2  C  p  ^2v  cx-*Lqw  2  e  ■* 

P^oo  oo 


(3.15) 


(3.16) 


(3.17) 


From  Fig. 3. 8  where  the  shear  stress  plot  as  a  function  of  t 
is  given,  the  heat  transfer  at  the  wall  can  also  be  obtained  from  the  re¬ 
lation 


02 

a-j  Tw 


(3.18) 


CHAPTER  IV 


VISCOUS  INTERACTIONS 

4.1  The  Displacement  Thickness 
Using  the  relation 


v-[peV  -  Jo“(Pe0e  -  pp)dy] 

+  ft  ^peA  '  Ope  "  p'dy-l  =  0  (4.1) 

for  the  unsteady  boundary  layer  as  given  in  reference  [19]  one  may  write, 
for  constant  pg  and  constant  Ug  , 

ft  <A-V  +  Ue  h  <A’6*>  ■  0  <4-2> 

where 

«*  =  /e(l  -  (4.3) 

0  peU 

e 

and 

6p  ■  /eO  -  ^)dy  (4.4) 

K  0  we 

From  expression  (3.8),  one  may  obtain 
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9-1  =  [9*  +  ^  "e2f  HI-'1  >  '  0-f) 


(4.5a) 


and 


g-f'  =  [gw  +  ^-Me2f'](i-f') 


(4.5b) 


Also  for  the  case  of  no  heat  transfer  g 1 ( 0 ,t )  =  0  and  from 
expression  (3.7),  one  may  write 


9  =  1 


+  M  2 

2  e 


(4.6) 


Using  expressions  (2.7),  (2.14),  (4.5a),  (4.5b)  and  (4.6)  we 
may  write,  following  a  procedure  outlined  in  Ref.  [3],  for  the  case  of 
no  heat  transfer  at  the  wall 


6*  = 


4^—]  J*  [1  +  ilTLi-  Me2  (l+f')](l-f')dn 


(4.7) 


and 


2v  cx 

/r  00 

6p  = 


•]  /e  {[1+  M  2(l+f')](l-f')  -  (l-f ‘ )}dn  (4.8) 

0 


Similarly  for  the  case  of  non-zero  heat  transfer  at  the  wall  we 


obtain 


2v  cx 

6*  -  4-g-]  /  [gw  * 

e  U 


i^U-Me2f'](l-f')dn 


(4.9) 
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and 


2v  cx  (  i\  o 

6p  =  4-g - ]  Je  {[gw  +  ilfii  MeZf'](l-f')  -  (l-f '  )}dn  (4.10) 

For  the  steady  state  case  dependence  on  time  would  disappear 
and  as  such  only  the  expressions  (4.7)  and  (4.9)  would  be  required  with 
the  function  f'  replaced  by  its  final  steady  state  value  f£  . 

4.2  Steady  State  Boundary  Layer  Thickness 

In  order  to  justify  the  Prandtl  Boundary  Layer  simplifications, 
calculations  were  made  for  6/x  for  the  initial  and  the  final  steady  state. 
By  definition 

p  U  (6-6*)  =  Je  p  u  dy  (4.11a) 

e  e  0 

or 

s-6*  =  /e  -frr  (4.iib) 

0  pee 


In  terms  of  the  transformed  variables 


6 

x 


^  2v  cx  a 

f  +  1  /[-n— ]  J6  f ' 

0 


dn 


(4.12) 
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This  result  may  also  be  expressed  in  terms  of  the  Hypersonic 
similarity  parameter  x  [see  equation  (4.19)]  as: 


(4.13) 


Figure  (4.1)*  shows  the  plot  of  6/x  as  a  function  of  x 
for  the  selected  cases  of  zero  and  non-zero  heat  transfer  at  the  wall. 


r  _ 

A  plot  of  ~/ x  as  a  function  of  M  (5<M  <20)  has  also  been  given  in 

A  c  c 


Fig.  (4.2)  for  the  case  of  zero  heat  transfer  at  the  wall. 

4.3  Weak  Interaction  Induced  Pressure 

In  the  acoustic  approximation,  the  pressure  rise  at  the  plate 
surface  and  the  outward  velocity  of  the  "displacement  plane"  may  be 
written  [17]  as,  respectively. 


(4.14) 


and 


(4.15) 


From  (4.2),  (4.14)  and  (4.15)  we  may  also  write 


*Curves  for  CHAPTER  IV  have  been  included  in  Appendix  ' F1. 
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Peae 


35 
'  3 1 


+  U 


35*  v 
e  3x  ' 


or 


=  Y(t 


35 

at 


M  M*) 
e  3x  ' 


Differentiation  of  expression  (4.9)  gives  [3] 


^  <2??  +  I?) 


3x 


(4.16) 


(4.17) 


where 


1  =  Jo[9w  +  ili^Me2f,](1-f')dn  (4-18) 

for  the  case  of  non-zero  wall  heat  transfer.  For  the  case  of  no  wall 
heat  transfer,  expression  (4.6)  should  be  used  in  (4.18) 

Using  the  hypersonic  similarity  parameter  x  [19]  defined  as 

c 

X  =  Me3/Cr?^M  e3/[Ri~]  <4-19> 

G  X  6 

and  noting  that  I  in  expression  (4.18)  is  a  function  of  x  only,  we 
get 
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M  14*  = 

e  3x 


/2  M. 


X 


(I-2t  |I) 


(4.20) 


Similarly  differentiation  of  expression  (4.10)  gives 


9J 

at 


(4.21) 


where 


J  =  Joe{[9w  +  He2f,](1-f')  '  (4.22) 

for  the  case  of  non-zero  wall  heat  transfer.  For  the  case  of  no  wall 
heat  transfer  one  should  use  expression  (4.6)  in  (4.22). 

Noting  expression  (4.19)  and  that  J  is  a  function  of  x 
only,  one  may  also  write 


/2  -  9J 

ae  3t  “  ^2  X  9t 


(4.23) 


Substitution  of  (4.20)  and  (4.23)  in  (4.16)  yields 


X 


[/2|i  + 


1 

72 


(I-2t  |I)]  X 


(4.24) 


Expression  (4.24)  has  been  derived  in  detail  in  Ref.  [3], 
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In  Ref.  [26],  parameter  X  has  been  used  in  place  of  x» 
where  X  and  x  are  related  by 


A  =  aX  (4.25) 

In  the  above  expression  "a"  depends  primarily  on  wall  temperature  for 
the  steady  state  case. 

Our  result  of  expression  (4.24)  may  also  be  expressed  in  terms 
of  the  parameter  X  as: 


Ap 

r-^=X  =  aX  (4.26) 

pe 

Here  "a"  depends  on  g  and  x  for  given  y  and  M  . 
Expression  (4.26)  is  similar  to  expression  (8)  of  Ref.  [27] 
where  the  steady  state  induced  pressure  results  have  been  obtained. 

For  the  final  steady  state  we  may  write,  from  expression  (4.24), 
as 


Apv2  yI2  - 

-T —  =  - ?  x 

pe  /2 

e 


(4.27) 


where 


Me2f'](l-f-)dn 


(4.28) 


- 


for  the  case  of  non-zero  wall  heat  transfer.  Replacement  of  gw  in 
(4.28)  by  expression  (4.6)  is  esseptial  for  the  case  of  no  wall  heat 
transfer. 

Finally,  one  obtains  the  transient  contribution  to  the  in¬ 
duced  pressure  referenced  to  the  final  steady  state  as: 

aj,  aPy2  -  flPy 
Pe  Pe 

=  [?l  •  /2  §7  ■  7T  (i_2t  §7)]  X  (4‘ 

e 

This  follows  directly  from  expressions  (4.24)  and  (4.27). 

Numerical  results  for  the  transient  contribution  to  the  pres 
sure  distribution  have  been  obtained  by  programming  expression  (4.29) 
on  IBM-360/67  Computer.  Various  results  are  displayed  in  Figs.  (4.3) 
through  (4.8)  in  Appendix  'F1. 


CHAPTER  V 


DISCUSSION  AND  CONCLUSIONS 


5.1  Discussion  of  Results 

(1)  For  zero  and  non-zero  wall  heat  transfer  transient  values  have 
been  obtained  for  the  shear  stress  at  the  wall  and  the  induced  pressure 
distribution  in  the  region  downstream  of  x  =  U  t/1.3  (0<t<1.3).  The 
selected  results  are  given  in  Appendix  'F1. 

(2)  The  stream  function  increment  ( Af )  distribution  has  been  given 
in  Figs.  3.3  and  3.4.  These  distributions  have  been  obtained  by  taking 
one  percent  change  in  the  free  stream  velocity.  The  discontinuity  in 

Af  -  distribution  is  noticeable  approximately  at  t  =  1.85  for  n  =  0.3, 
t  =  1.65  for  n  =  1 ,  x  =  1.4  for  n  =  2,  x  =  1.3  for  n  =  5. 

It  should  be  mentioned  that  Af  -  distribution,  for  the  stepwise- 
accelerated  and  decelerated  flat  plate  cases,  happens  to  be  identical. 
Therefore,  the  results  of  Ref.  [3]  may  be  extended  to  the  range  1<t<1.3. 

(3)  The  highest  order  derivative  terms  (in  time  and  space)  of 
equation  (2.31)  are 


'  '  2(l-xf£)  ^  '  +  — 


(5.1a) 


or 


3(Af 1 ) 

8t 


K 


32(Af] 

3n2 


+  — 


(5.1b) 
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Expression  (5.1)  has  the  character  of  the  conduction  heat  equation  with  the 
coefficient  of  conductivity  given  by 

K  =  erPFfp"  •  (5'2) 

In  expression  (5.2)  K  would  have  different  signs  for  0<T<ir  ,  and 

1  1  c 

it  <t<°°  .  At  t  =  7T  »  K  would  be  infinite.  It  appears  that  the  dis- 

t2  t2 

continuity  in  the  Af  -  distribution  occurs  upon  approaching  the  neighbour¬ 
hood  of  t  =  ir  .  A  line  denoting  t  =  ir  is  shown  in  Fig.  (3.3). 
t2  t2 

(4)  Some  attempts  were  made  (see  Appendices  'A',  ' B '  and  1 D 1 )  to 
use  the  11-point  molecule  in  place  of  the  7-point  molecule.  The  11-point 
molecule  is  the  result  of  the  commonly  used  expressions  for  central  de¬ 
rivatives.  The  results  are  shown  in  Fig.  A. 2  and  Fig.  B.2.  It  has  been 
found  that  the  methods  incorporating  the  11-point  molecule  are  not  stable 
for  the  equation  (2.31)  with  the  associated  time  and  boundary  conditions 
given  in  section  2.5.  The  Af  -  distribution  of  Fig.  3.3  and  Fig.  3.4  is 
based  on  the  recurrence  relations  (3.1),  (3.2),  and  (3.3)  incorporating 
the  7-point  molecule.  Double  precision  was  used  in  the  computations. 

Fig.  3.4  shows  that  Af  is  a  monotonical ly  increasing  function  of  n 
for  fixed  value  of  t  . 

(5)  Figures  3.5,  3.6  and  3.7  indicate  the  satisfaction  of  the  dif¬ 
ferential  equation  (2.31)  by  the  stream  function  increment  Af  and  its 
derivatives.  The  quantity  1 d 1  defined  by 

d  =  2(Af ) '  +  2T[f“(Af)  -  f  2  ( Af )  1  ]  -  (Af)"' 


-  f2(Af)"  -  f£'(Af) 


(5.3) 
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has  been  calculated.  The  Af  -  distribution  has  been  obtained  using  the 
7-point  molecule.  But  in  the  computation  of  the  quantity  ' d 1 ,  the  de¬ 
rivatives  based  on  the  11-point  molecule  are  used.  The  discrepancies 
(d^O)  are  the  result  of  this  change-over  from  7-point  molecule  to  the 
11 -point  molecule.  The  satisfaction  of  the  differential  equation  is 
quite  reasonable  for  large  values  of  n  and  t  .  The  discrepancies 
(d^O)  are  comparable  with  the  highest  order  derivative  ( Af ) ' ' '  for 
smaller  values  of  r\  and  t. 

As  we  mentioned  previously,  the  numerical  methods  incorporating 
the  11-point  molecule  have  been  found  unstable  for  the  present  problem. 

It  is  felt  that  the  results  obtained  by  using  the  7-point  molecule  are 
quite  satisfactory  for  the  estimation  of  the  transient  shearing  stress 
and  the  induced  pressures. 

(6)  In  Fig.  3.8  the  transient  distribution  of  Af''(0,T)  is  shown. 
This  gives  immediately  the  desired  shear  stress  distribution  at  the  wall. 
For  large  values  of  t  ,  Af 1 1 ( 0 ,t)  appears  to  approach  the  final  steady 
state  value. 

We  noted  in  section  3.4  that  the  heat  transfer  at  the  wall  is 
related  to  the  shear  stress  at  the  wall  through  the  Reynolds  Analogy: 


0*2 

a-|  Tw 


(5.4) 


Therefore,  heat  transfer  at  the  wall  may  also  be  obtained  from  Fig.  3.8. 

(7)  The  transient  contribution  to  the  induced  pressure  distribution 
for  the  case  of  no  heat  transfer  at  the  wall  is  given  in  Fig.  4.3a.  The 


.  ■ 
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transient  contribution  becomes  laroer  at  smaller  Mach  Numbers  and  at 
bigger  values  of  t  .  The  absolute  magnitude  of  the  transient  contri¬ 
bution  correspondi na  to  Mg  =  5  is  almost  equal  at  t  =  1.0  and  larger 
at  t  =  0.1  as  compared  with  the  correspondi na  contribution  for  the  step¬ 
wise-accelerated  flat  plate  of  Ref.  [3].  This  may  be  noticed  from  the  com¬ 
parison  of  Figs.  4.3a  and  4.3b.  Also  the  contribution  is  larger  at 
Mg  =  30  and  x  =  0.1  and  x  =  0.05  in  this  case. 

In  Fiq.  4.3c  /x  has  been  plotted  as  a  function  of  x. 
pe 

This  is  a  precise  form  of  Fig.  4.3a.  One  may  notice  abrupt  convexities 
after  x  =  1.0  in  this  figure.  These  convexities  may  indicate  that  we 
are  soon  approachinci  the  vicinity  of  the  "blow-up"  reaion  of  Af  -  distri¬ 
bution. 


(8)  Fiqures  4.4a,  4.5a  and  4.6a  qive  the  transient  contribution 
to  the  pressure  distribution  with  non-zero  heat  transfer.  In  Fiqures  4.7 
and  4.8,  some  of  these  results  have  been  further  analysed  for  x  =  3.5, 

Mg  =  10  and  20  respectively.  It  is  noticed  that  at  a  fixed  Mach  Number 
and  fixed  x  the  transient  contribution  is  larger  for  lower  values  of  q 

w 

up  to  x-0.18  whereas  the  transient  contribution  is  smaller  for  lower 

values  of  q,  in  the  ranqe  0 . 1 8<x< 1 . 3 .  For  the  case  of  an  adiabatic  wall 
w 

the  transient  contribution  to  the  pressure  distribution  is  nearly  constant 
as  a  function  of  x  up  to  x  =  1;  after  x>l  the  transient  contribution 
increases . 


The  absolute  magnitude  of  these  contributions  is  found  to  be 
larger  in  the  range  0<_gw<l  and  smaller  for  gw>_l  as  compared  to  the  corre¬ 
sponding  transient  contributions  for  a  stepwise-accelerated  flat  plate. 
Fiqures  4.4b,  4.5b  and  4.6b  of  Ref.  [3]  have  been  included  for  the  sake 


' 

■  x 


43 


of  comparison. 

(9)  For  a  given  g  ,  given  x  anc*  given  t  ,  the  transient  con- 

w 

tribution  decreases  with  the  increase  in  Mach  Number.  This  may  be  noticed 
again  from  Figures  4.4a,  4.5a  and  4.6a. 

For  large  Mach  Numbers  and  large  values  of  t  the  transient 
contribution  becomes  independent  of  the  Mach  Number  as  gw  approaches 
zero.  Figure  4.7  and  4.8  along  with  Figs.  4.4a,  4.5a  and  4.6a  may  be 
analysed  to  indicate  this. 

Figures  4.4c,  4.5c  and  4.6c  have  been  included  to  reduce  the 
number  of  curves  of  the  respective  Figs.  4.4a,  4.5a  and  4.6a. 

(10)  Now,  consider  the  expression  (4.29): 

=  ff?  [?2  (I2_I)  '  /2  (l?  '  1  §7)]  *  (5'5) 

e 

For  the  case  of  no  heat  transfer  at  the  wall  and  0<r <_1 . 3 ,  it 
has  been  found  that 


1 

72 


a2-n 


>> 


72  ( 


9J 

9t 
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and,  therefore,  one  may  use  the  approximate  relation 


M  = 
Pe 


_ 1 _ 

72  M  2 
e 


(i2-d 


X 


(5.6) 


which  is  consistent  with  the  findings  of  Ref.  [3], 
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For  l<g,  <81  the  contribution  of  all  terms  must  be  considered. 
— Jw— 

A  I 

However,  for  g  <1  and  t«1  ,  it  has  been  noticed  that  /2  is  the 

W  oT 

predominant  term  and  for  this  case  we  may  use  the  approximate  relation 


Ap  /2  y  9J  — 
p_  3t  x 


(5.7) 


This  result  is,  also,  consistent  with  the  corresponding  stepwise-accelerated 
flat  plate  result  of  Ref.  [3].  The  solution  in  the  region  t-K),  should 
correspond  to  that  given  by  Lord  Rayleigh  [8]. 

(11)  Finally,  we  justify  the  use  of  the  Prandtl's  Boundary  Layer 
simplifications  in  the  present  work.  As  pointed  out  by  Shen  [1],  one 
may  use  the  Prandtl's  Boundary  Layer  equations  for  the  Hypersonic  Flows 

r  o  r 

provided  (7)  «1.  Calculations  of  7  have  been  made  for  the  initial  and 

a  A 

the  final  steady  state  cases.  Some  of  the  limiting  results  are  presented 

6  2 

in  Figs.  4.1  and  4.2.  It  has  been  found  that  the  quantity  (7;  never  ex- 
ceeds  0.3  in  the  weak  interaction  rage  (i.e.  x  =  0  to  X  -  3.5). 


5.2  Concluding  Remarks 

The  results,  presented  in  the  previous  section,  have  been  ob¬ 
tained  for  the  initial  zero  pressure  gradient  case  and,  therefore,  are 
valid  for  weak  interaction  induced  pressures.  The  numerical  results 
cover  the  range  0<r<_  1 . 3  and  additional  attempts  should  be  made  to  obtain 
the  solution  in  the  region  1<t<?°.  It  would  be  convenient  to  use  a  para¬ 
meter  a  defined  by 
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(5.8) 


in  this  range.  The  solution  at  t=of1  may  serve  as  the  additional 
boundary  condition  for  the  region  in  question. 

Once  the  main  aspects  of  the  unsteady  flow  over  a  flat  plate 
at  zero  angle  of  attack  are  known  for  0<r<_°°  ,  it  would  seem  interesting 
to  extend  the  analysis  to  include  the  cases  of  a  wedge,  an  inclined  flat 
plate,  and  the  circular  cone.  Further,  the  problems  of  strong  interaction 
induced  pressure  and  the  oscillating  airfoils  may  be  treated. 
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APPENDIX  'A 


RECURRENCE  RELATIONS  INVOLVING  THE  11 -POINT  MOLECULE 

A. 1  Formation  of  the  Difference  Equation 

For  zero  pressure  gradient  we  have  from  equation  (2.31): 

2(Af)'  +  2x[f“(Af)  f  ^  ( Af ) ']  -  Af'" 

-  f2(Af)'  '  -  "f 2 *  ( Af )  =  0  (A.  1 ) 


or 
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(A. 2) 
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Let  H  and  T  be  the  elementary  step  magnitudes  in  n  and  x 
directions  respecti vely .  Then  [28]  we  may  write 


^  =  2T  (Sm+1  ,n  "  Sm-1  ,n^ 


(A. 3) 


S'  '  2H  (Sm,n+1  '  Sm,n-1* 


(A. 4) 


S'1  =  -U-  (S  ,  -  2S  +  S  ) 

1^1^-  m,n+ 1  m$n  m ^ n “  I 


(A. 5) 


2H 


T  (Sm  n  i  o  -  2S  , -i  +  2S  -t  -  S  9) 
3  m,n+2  m,n+ I  m,n-l  m,n-2 


(A. 6) 


S'  4HT  (Sm+l,n+l  "  Sm-1 ,n+l  +  Sm-1 ,n-l  "  Sm+l,n-l^ 


Substituting  (A. 3)  through  (A. 7)  into  expression  (A. 2),  we 


obtai n 


2 

H  (Sm+l,n+l  "  Sm-l,n+l  +  Vl  ,n-l  "  Sm+l,n-l^ 


+  4  Tm,n  H3T^f2m,n  2T  ^Sm+l,n  '  Sm-l,n)  "  f2m,n  4HT  ^Sm+1 ,n+l 


"  Sm-l,n+l  +  Sm-1  ,n-l  "  Sm+l,n-l^  '  T^Sm,n+2  "  2Sm, 


n+1 


+  2Sm,n-l  "  Sm,n-2^  "  ^2m,n  2HT^Sm,n+l  “  2Sm,n  +  Sm,n-1^ 


-  f'  '  n  2H  T  Sm  =  0 
2m, n  m,n 


(A. 8} 
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A.  2  Derivation  of  the  Recurrence  Relations 


form 


where 


Expression  (A. 8)  can  also  be  written  in  the  following  general 


am,n  ^m+1 ,n+l  +  ^m,n  ^m+1 ,n  +  cm,n  ^m+1 ,n-l  ^m.n  (A. 9) 


a  =  (1-t  f9  )  H‘ 
m,n  m,n  2m, n 


(A. 10) 


b  n  =  2f '  '  x  H' 
m,n  2m,  n  m,n 


(A. 11) 


C  =  “  (1-T  f«  )  H‘ 
m,n  m,n  2m, n 


(A. 12) 


d  =  T(S  ,9  -  2S  +  2S  -  S  9) 
m,n  m,n+2  m,n+l  m,n-l  m,n-2 


+  2HT f  n  n  (Sm  , -i  -  2Sm  +  Sm  n  ,) 
2m,  n  m,n+l  m,n  m,n-l 


+  Vn  +  <Vl,n+l  ’ 


+  2H3t  f '  '  S  ,  „ 
m,n  2m, n  m- I ,n 


(A. 13) 


Considering  the  boundary  conditions  in  conjunction  with  (A. 9)  a  recurrence 
formula  can  be  obtained,  namely, 
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lm,n  ^m+1  ,n+l  +  ^m,!! 


S  =  d 1 
m+1 ,n  m,n 


(A. 14) 


where 


k,  ,  m,n-l  ~ 

b  =  b  -  c  r-r2 -  ;  n>3 

m,n  m,n  m,n  b^  _1 


(A. 15) 


and 


d'  =  d  -  c  JMld.  . 
m,n  m,n  m,n  b'  , 

II I  5  I  I  ""  I 


n>3 


(A. 16) 


and  where  for  n  =  3  we  have  from  the  boundary  conditions 


b 1  =  b  o+ic  o 

m,3  m,3  4  m,3 


(A. 17) 


d '  0  =  d  0 
m,3  m,3 


(A. 18) 


Having  initial  b'  „  ,  and  d'  ,  we  can  find  from  (A. 15)  and 
3  m,n-l  m,n-l 

(A. 16)  all  b^  n  and  d^  n  respectively .  With  these  values  we  can  obtain 
from  (A. 14)  variable  Sm+-|  n  provided  that  we  have  Sm+-|  n+-j  , 

Let  n  =  N  represent  the  location  for  which  "e"  condition  pre¬ 
vails.  We  know  at  this  location  from  the  north  boundary  condition  that 


3S 

3n 


2ST  <-3Si 


+  4S2  -  S3)  =  0 


(A. 19) 
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Therefore, 


"  3Sm+l  ,N  +  4Sm+l ,N-1  '  Sm+1  ,N-2  ‘  0 


(A. 20) 


Also  from  (A. 14)  one  may  write 


am,N-l  Sm+1,N  +  bm,N-l  Sm+1  ,N-1  "  dm,N-l 


(A. 21a) 


and 


am,N-2  Sm+1,N-1  +  bm,N-2  Sm+1 ,N-2  “  dm,N-2 


(A. 22a) 


From  (A. 22a)  one  obtains 


d‘  a 

-  _  m,N-2  dm,N-2  - 

m+1,N-2  bm,N-2  bm,N-2  m+1 ’N“1 


(A. 22b) 


Substituting  (A. 22b)  into  (A. 20)  we  have 


or  _  ac  +  m,N-2  _  m,N-2  _  n 

m+l,N  m+l,N-l  b^N_2  b^N_2  m+1  ,N-1 


(A. 23) 


From  (A. 21a)  we  have 


m,N-l  _  m,N-1 

5m+l  ,N-1  b '  m  -i  "  b '  m+1  ,N 

m,N-l  m,n- 1 


(A. 21b) 
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Substituting  (A. 21b)  into  (A. 23)  one  finally  obtains 


+ 


(A. 24a) 


or 


bm,N-l  bm,N-2  bm,N-2  bm,N-l 


S 


m+l  ,N 


3  +  4  Vn-1  +  am,N-2  am,N-l 


(A. 24b) 


bm,N-l  bm,N-2  bm,N-l 


The  above  expression  constitutes  the  starting  point  in  the  use  of  ex¬ 
pression  (A. 14). 


Expressions  (A. 14)  through  (A. 18)  and  (A. 24)  were  programmed  on 


IBM-360/67  computer.  Results  are  shown  in  Fig.  A. 2.  Values  of  Af 
for  t  =  0  have  been  calculated  from  the  initial  condition  (2.37)  pre¬ 
scribed  for  the  problem.  Recurrence  relations  (3.1),  (3.2),  and  (3.3) 
have  been  used  to  calculate  Af  for  t  =  0.05  (or  m=2).  For  x  =  0.1 
(or  m=3)  Af  has  been  calculated  using  the  recurrence  relations  (A. 14) 
through  (A. 18)  and  (A. 24).  After  t  =  0.1,  the  method  described  in  this 
appendix  becomes  unstable.  The  values  of  Af  ,  obtained  for  x  =  0.1, 
are  about  0.15%  higher  (for  n  =  I'.O  and  n  =  5.0)  than  the  corresponding 
values  obtained  by  using  the  recurrence  relations  (3.1),  (3.2),  and  (3.3). 


55 


Fig.  A. 2  Stream  Function  Increment  (Af)  Distribution  Curve 
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A. 3  Reasons  for  the  Unstability 

When  a  search  was  made  to  find  out  the  causes  for  the  unstability 
it  was  noticed  to  be  associated  with  the  expressions  (A. 15)  through  (A. 18). 
From  expression  (A. 17)  it  may  be  noted  that  the  coefficient  of  cm  g  is  1/4 
Also  one  may  write  from  (A. 15)  and  (A. 17)  as: 


b  n  b  /]  "I"  C  n  ( — 

m,4  m,4  m,4 


1  o 

m,3 


m,3  4  m,3 


) 


(A. 25) 


In  the  above  expression  b  and  b  „  are  of  the  order  of  0.000095 

r  m,3  m,4 

or  smaller  for  t  =  0.10  whereas  am  Q  and  c  ,  are  of  the  order  of 

m,J  m,  J 

0.00995  or  smaller  and  c  ,  =  -  a  0  .  For  a  rough  estimate  of  the 

m,3  m,3  3 

coefficient  of  c  n  ,  we  may  drop  bm  Q  and  we  find  the  coefficient  .. 

m,4  r  m,3 

to  be  4. 


Again,  one  may  obtain  from  expressions  (A. 15)  and  (A. 25) 


b'r=b  r+C  r  ( - 

m,5  m,5  m,5 


am,4*3m,3  +  4  am,4cm,3 


b  nb  .  +  -;<•  b  n  c  —  a  n 
m,3  m,4  4  m,4  m,3  m,3  m,4 


-)  (A. 26) 


Here  also  bm  ,  bm  „ 
m,3  m,4 

c  ,  ,  c  „  ,  am  o  and 
m,3  m  ,4  rn  $  3 


and  b  r  are  of  the  order  of  0.000095  whereas 
m,5 

a  „  are  of  the  order  of  0.00995  with  c  =-am  „ 
m,4  m,n  m,n 


Neglecting  the  smaller  order  terms  we  may  obtain  the  coefficient  of  c  r, 

m )  o 

approximately,  as 


1  am,4  cm,3 
4  am,3  cm,4 


57 


which  gives  1/4. 

Thus,  we  find  that  the  coefficient  of  c  „  in  the  recurrence 

m,n 

relation  (A. 15)  fluctuates  between  1/4  and  4.  This  causes  fluctuations 

in  the  value  of  b'  „  and  d1  M  and  consequently  the  S  „  values  obtained 

m,n  m,n  ^  m,n 

from  (A. 14)  fluctuate.  There  are  no  fluctuations  in  the  values  of  a 
It  was  felt  that  the  unstability  in  using  the  full  molecule 
may  be  associated  with  the  recurrence  relations  used  and,  therefore,  the 
Gauss's  Elimination  method  was  tried.  This  has  been  described  in  Appendix 
' B ' .  However,  the  unstability  was  noticed  with  the  Gauss's  Elimination 
method  also.  A  third  method  described  in  Appendix  ' D '  did  not  give  any 
satisfactory  results  either.  It  is  thought  that  the  11-point  molecule 
is  not  suitable  for  the  computations  of  Sm  n  of  equation  (A. 9). 
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APPENDIX  ' B ' 


GAUSS'S  ELIMINATION  METHOD 


B.l  Formation  of  the  Problem  in  a  Tri diagonal  System 


t 
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Fig.  B.l 


The  Gauss's  Elimination  Method  has  been  used  to  solve  the 


equation 


\n,n  $m+l  ,n+l  + 


b  S  +  c 
m,n  m+1 ,n  m,n 


^m+l,n-l  "  dm,n 


(B.l) 


derived  in  Appendix  'A'  [see  equation  (A. 9)]. 

The  associated  boundary  conditions  are  (see  section  2.5): 
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as 


Ffl+  1  j  1  _  n  _  r 

an  m+l,l 


and 


(B.2) 


as 


m+l  ,100  n 

3n 


In  addition  we  have  the  initial  condition  prescribed  by  (2.37). 

The  first  of  the  two  south  boundary  conditions  implies 


2An  3Sm+l  ,1  +  4Sm+l  ,2  "  Sm+1 ,3^  "  0 


and  using  the  second  south  boundary,  Sm+-|  -j  =  0,  we  have 


Sm+1 ,2  =  4  Sm+1 ,3 


(B.3) 


With  n  =  3,  we  obtain  from  equation  (B.l) 


am,3  Sm+1 ,4  +  bm,3  Sm+1 ,3  +  cm,3  Sm+1 ,2  =  dm,3 


In  view  of  relation  (B.3),  it  becomes 


a  o  S  ,-,  „  +  (b  0+tc  ?)  S  ^  =  d  -3 
m,3  m+l  ,4  v  m,3  4  m,3'  m+l  ,3  m,3 


(B.4) 
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With  n  =  4,5,0..,  97,  98,  one  obtains  from  (B.l) 

am,4  Sm+1 ,5  +  bm,4  Sm+1 ,4  +  cm,4  Sm+1 ,3  =  dm,4 

am,5  Sm+1 ,6  +  bm,5  Sm+1 ,5  +  cm,5  Sm+1 ,4  =  dm,5 

UB.5) 

am,97  Sm+1 ,98  +  bm,97  Sm+1 ,97  +  cm,97  Sm+1 ,96  =  dm,97 

am ,98  Sm+1 ,99  +  bm,98  Sm+1 ,98  +  cm,98  Sm+1 ,97  =  dm,98  J 

The  north  boundary  condition  gives 

Sm+1 ,100  =  I  ^4Sm+l  ,99  "  Sm+1 ,98^ 
and  expression  (B.l)  with  n  =  99  gives 

am,99  Sm+1 ,100  +  bm,99  Sm+1 ,99  +  cm,99  Sm+1 ,98  =  dm,99 
Using  (B.6) ,  we  obtain 

am,99  +  bm,99^  Sm+1 ,99  +  (cm,99  "  3  am,99^  Sm+1 ,98  =  dm,99 
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In  this  manner  we  obtain  97  equations  and  97  unknowns,  namely, 

$m+l  3**’**  $m+l  99’  Putt"*n9  the  equations  (B.4),  (B.5)  and  (B.7)  in 
the  matrix  form,  we  find 


—  — 

bm,3  +  4  cm,3  am,3 

Bm+1 ,3 

dm,3 

cm,4  bm,4 

am,4 

Sm+1 ,4 

dm,4 

Sm+1 ,5 

• 

• 

• 

• 

• 

= 

• 

• 

m  < 

• 

Sm+1 ,97 

dm,97 

Cm,98  bm,98 

am,98 

Sm+1 ,98 

dm,98 

c  -  1  a 

m,99  3  m,99 

4 

3  am,99  +  bm,99 

Sm+1 ,99 

dm,99 

B.2  Gauss  Elimination  for  a  Tri diagonal  System 

In  the  previous  section  we  obtained  a  system  of  97  linear  al¬ 
gebraic  equations  of  the  form 


B3  ^3  ^4  =  ^3 


Ai  Si-1  +  BiSi  +  Ci  Si+1  =  Di ; 


(B.9) 


Am  Sm  t  +  BmSm  =  Dm;  N  =  99 
N  N-l  N  N  N 


(B.8) 


62 


where  the  S.  are  unknown  and  where  the  A.  ,  B.  ,  C.  and  D.  are  given 
coefficients.  Since  the  coefficient  matrix  of  the  system  is  tridiagonal, 
i.e.  has  no  non-zero  elements  except  on  the  main  diagonal  and  on  the 
diagonals  adjacent  to  the  main  diagonal,  there  is  a  very  simple  algorithm 
for  solving  the  equations.  It  was  first  used  in  connection  with  para¬ 
bolic  partial  differential  equations  by  Bruce,  Peaceman,  Rachford  and 
Rice  [24],  In  this  algorithm,  which  is  actually  just  the  Gauss's  Elimi¬ 
nation  method,  the  number  of  arithmatic  operations  is  proportional  to  N 

1  3 

as  compared  to  y  N  for  the  general  matrix. 

The  S .  as  defined  by  expression  (B.9)  can  be  written  out  ex¬ 
plicitly  through  use  of  the  following  formulas: 


^3  ^i 

b*  =  /  ,  b*  =  ft  ■  a1-  ■■  ,  i=4,5, . . .  ,N-1 

J  b3  1  i  i  i-1 


D3  ^i"^i  ^i-l 

0*  =  —  n*  =  — - ! - ! — -  i=4  5  N-l 

q3  B3  »  qi  B.-Ai  b*_-,  *  1  H,D,,0”IN  1 


(B. 10) 


With  these  coefficients  the  unknown  S.  are  given  by 


=  qN 


(B.ll) 


Si  =  qi"bi  Si+1  »  1  N_1  »N-2,-  •  •  >3' 


The  values  of  S.  obtained  by  using  this  particular  method  are 
shown  in  Fig.  B.2.  Initial  condition  (2.37)  prescribes  the  value  of 
S.(or  Af)  for  t=0.  For  t=0.05  and  t=0.1,  the  value  of  S.  have  been 
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Fig.  B . 2  Stream  Function  Increment  (Af)  Distribution  Curve 
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obtained  from  the  recurrence  relations  (3.1),  (3.2),  and  (3.3).  Gauss's 
method  gives  the  value  for  t>0„15.  We  notice  from  the  figure  that  the 
method  of  this  appendix  becomes  unstable  after  t=0.15. 


APPENDIX  'C 


ITERATIVE  FINITE  DIFFERENCE  METHOD 
INVOLVING  THE  7-POINT  MOLECULE 


C. 1  Formation  of  the  Difference  Equation 

The  iterative  finite  difference  method  has  been  developed  to 
solve  equation  (2.31).  The  finite  difference  approximations  involve  the 
7-point  molecule. 


T - i 

P — 

>N 

/  V  .  .X 

,N“1 

v  r  ■ j 

N-2 

>n+l 

•  \  .  * 

n-1 

— < 

pi 

m  -1  m  m  ♦  1 


t  — — 


Fig.  C.l 


If  H  and  T  denote  elementary  step  magnitudes  in  n  and  x 
direction,  then  one  may  write  [29] 
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=  j  [Af(m,n)  -  Af(m-1  ,n)]  (C.l) 

f~-  =  hi  tAf(m»n+1)  ‘  Af (m,n-l ) ]  (C.2) 

2 

|^~  =  [Af(m,n+1 )  -  Af(m,n-1)  -Af(m-l,n+l) 

+  Af (m-1 ,n-l )]  (C.3) 

2 

— -a—  =  -U-  [Af(m,n+1)  -  2Af(m,n)  +  Af(m,n-1)]  (C.4) 

3n  n 

and  using  Taylor's  Series 

3 

=  -U-  [Af(m,n+1)  -  3Af(m,n)  +  3Af(m,n-l)  -  Af(m,n-2)]  (C.5) 

3n  r 

Substitution  of  expressions  (C.l)  through  (Co5)  into  (2.31) 

yields : 

T(Af(m,n+l)  -  3Af(m,n)  +  3Af(m,n-l)  -  Af(m,n-2)]  +  HTf2(m,n)  X 

[Af(m,n+1)  -  2Af(m,n)  +  Af(m,n-1)]  +  H3Tf2 1 (m,n) Af (m,n) 

=  H2[l-Tmf2(m,n)]  x  [Af(m,n+1)  -  Af(m,n-1)  -  Af(m-l,n+l) 

+  Af (m-1 9 n- 1 )]  +  2H3Tmf2' (m,n) [Af (m,n)  -  Af(m-l,n)] 


(C.6) 
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Writing  (C.6)  in  the  general  form,  one  obtains 


(C.7a) 


or 


e  -  (a  S  , ,  +  c  S  i+d  S  „) 
m,n  v  m,n  m,n+1  m,n  m,n-1  m,n  m,n-2' 


m,n 


(C.7b) 


where 


am,n  =  T  +  HTf2(m’n)  '  H2[l-Tmf^(m,n)]  (C.8) 

bm,n  =  -  3T  -  2HTf2(m,n)  +  H3Tf“(m,n) 

-  2H3Tmf£'(m,n)  (C.9) 

cm,n  =  3T  +  HTf2(m,n)  +  H2[l-xmf^(m,n)]  (C.10) 


-  T 


em,n  '  H  ^1'Tmf2^n,’n*^Vl  ,n-l  "  Vl.n+l^ 


=  2H  Tmf-(m>n)  Sm.1 


(c.n) 
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The  associated  boundary  conditions  with  the  equation  (2.31)  are: 


(C. 12) 


(C . 1 3) 


S 


m,l 


=  0 


(C. 14) 


The  initial  values  of  S  „  may  be  obtained  from  (2.37). 

m,n  J 


In  the  above  expressions,  n  =  100  represents  the  location  for 
which  "e"  condition  prevails. 

C.2  Iterative  Finite  Difference  Method 

In  this  method,  one  assumes  the  values  of  the  function  and  then 
corrects  those  values  with  the  help  of  the  differential  equation  in  the 
finite  difference  form  and  the  associated  boundary  conditions. 


To  start  with,  one  may  assume  the  same  values  of  S  in  the 

m  ,n 


5 


second  column  of  Fig.  C.l  as  those  in  the  previous  column  [The  values  of 

S  in  the  first  column  will  be  given  by  the  initial  condition  ?2.37)]. 
m,n  3  J 


Using  (C.7b),  (C.12),  (C.13)  and  (C.14)  the  values  of  S  so  assumed  can 

III  )  1 1 

be  suitably  iterated  to  obtain  the  "exact"  values  of  the  function  S  in 

m  n 


m,n 


the  second  column  correct  to  the  desired  degree  of  accuracy.  In  obtaining 
the  values  of  Fig.  C.2,  error  was  driven  out  of  the  sixth  place  after  the 
decimal  point  in  the  value  of  S  .  The  computation  time  increases  rapidly 
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Fig.  C . 2  Stream  Function  Increment  (Af)  Distribution  Curve 
(Iterative  Finite-Difference  Method  of  Appendix  ' C ’ ) 
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as  we  want  to  drive  the  error  out  of  the  seventh  place  after  the  decimal 
point. 

Having  obtained  the  values  in  the  second  column,  one  may  proceed 
to  obtain  the  values  of  S  in  the  third  column.  This  marching  forward 
procedure  in  association  with  the  iterative  technique  may  be  used  for  the 
successive  columns  and  the  error  may  be  driven  out  of  all  the  8-places 
after  the  decimal  point  (in  single  precision)  depending  upon  the  accuracy 
desired. 

The  values  of  S  obtained  (see  Fig.  C.2)  using  this  method  are 

m,n  s/a 

about  4%  higher  (at  t=1 )  as  compared  with  the  values  of  Fig.  3.3*  The 
present  method  happens  to  be  costly  computer  time-wise. 

Finally,  it  should  be  mentioned  that  this  numerical  method  would 
prove  most  effective  for  functions  which  decay  more  rapidly  with  time  than 
the  function  under  consideration. 


*See  Appendix  'F 


APPENDIX  ' D ' 

ITERATIVE  FINITE  DIFFERENCE  METHOD  INVOLVING 
THE  11-POINT  MOLECULE 


An  iterative  finite  difference  method  similar  to  Appendix  'C' 
was  also  tried  for  the  full  molecule  consisting  of  11-points. 

Rewriting  equation  (A. 8)  of  Appendix  'A'  in  the  general  form  as 
shown  below,  one  finds 


A  S  ,  0  +  B  S  , +  C  S  +  D  S  -] 
m,n  m,n+2  m,n  m,n+l  m,n  m,n  m,n  m,n-l 


+  E  S  0  =  F 
m,n  m,n-2  m,n 


(0.1) 


where 


A  =  T 
m,n 


B  =  -  2T  +  2HT  f9m 
m  3  n  2m  3  n 


Cm  n  =  -  4HT  f9m  n  +  2H  T  f  “  n 
m,n  2m, n  2m, n 


D  =  2T  +  2HT  f0 
m,n  2m, n 


(D.2) 

(0.3) 

(D.4) 

(0.5) 
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Equation 


S 


E  =  -T 
m.n 


(D.6) 


^m,n  ^"Tm,n^2m,n^  ^m+1  ,n+l"^m+l  ,n-l  ^ 


+  2H3t  „  -  (1-t  f'  J  H2 

m,n2  m+ 1  ,n  m,n  2m,n 


x  (Sm  i  n+i  '  sm  l  „  i)  -  f“  S  ,  n  (D.7) 

m- I  ,n+ 1  m- I , n- I  m,n  2m, n  m-l,n 


(D.l)  may  also  be  written  as 


m.n 


^m,n  ^Dm,n^m,n-1 ~^m,n^m,n-2  ^m,n^m,n+2~^m,n^m,n+1 


m.n 


(D.8) 


The  associated  time  and  boundary  conditions  are  (see  Sec,  2,5) 


U 


U 


el 


’1  ,n>l 


f2 

e2  ^1 ,  U  ,  e2 

/[n^]n 
ue2 


+  (^  ~  D(n-1)H 


+  f. 


■l,n 


(D.9) 


_ m»100  =  o  =  — —  r-3  s  +  4  s  -  s  i 

3n  u  2An  L  m,100  *  in, 99  5m,98J 


(D. 10) 


8S  l  l 

— EL*  1  =  n  =  —  —  r_3  S  +  4  S  -  S  1 
9n  u  2An  L  J  m,l  ^  5m,2  \,3J 


(0.11) 
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and 


(D. 12) 


First,  the  S  distribution  was  obtained  for  0<t<1  using  the 
recurrence  relations  (3.1),  (3.2),  and  (3.3).  Later,  equations  (D.8) 
through  (D.12)  were  used  to  modify  the  Sm  n  distribution  to  incorporate 
the  full  molecule  into  the  calculations.  However,  the  convergence  (if 
any)  was  too  slow  and  the  method  was  given  up. 


APPENDIX  'E 


EFFECT  OF  THE  VARIATION  OF  GRID  SIZE 

Figure  E.l  shows  the  effect  of  the  variation  of  grid  size  upon 
the  value  of  stream  function  increment  Af  .  The  elementary  step  sizes 
in  the  n  and  t  directions  are  denoted  by  H  and  T  respectively.  As  we 
keep  the  value  of  H  fixed  and  vary  the  step  size  T  in  the  direction  of  x  , 
we  notice  that  the  maximum  variation  in  the  value  of  Af  happens  to  be 
about  5.8%  when  the  step  size  T  is  reduced  to  one-tenth  of  its  previous 
value  of  0.05.  However,  the  leading  derivative  terms  in  expression  (2.31) 
indicate  its  parabolic  nature  and,  therefore,  it  will  not  be  proper  to 
take  excessively  small  steps  in  the  T-direction,  which  will  lead  to  large 
round  off  errors  [25]. 

For  a  constant  value  of  T,  the  maximum  variation  in  the  value 
of  Af  with  the  variation  in  H  has  been  noticed  to  be  2.5 %  when  H  is  re¬ 
duced  to  one-half  of  its  previous  value  of  0.1. 

Lastly,  it  should  be  mentioned  that  the  smaller  step  sizes  also 
give  rise  to  computer-storage  problems  (especially,  if  one  wants  to  use 
double  precision  for  calculations)  and  also  the  numerical  approach  becomes 
costly  computer  time-wise.  Keeping  in  view  these  facts  it  was  decided  to 
keep  step  size  H  =  0.1  in  the  redirection  and  step  size  T  =  0.05  in  the 
x-di recti  on . 
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1/T 


Fig.  E.l  Effect  of  the  Variation  of  Grid  Size 
on  the  Value  of  Stream  Function  Increment  Af 


APPENDIX  ' F ' 

CURVES  FOR  CHAPTER  III  &  IV 
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Olx(jV)  cOlx(jV) 
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Fig.  3.3  Stream  Function  Increment  (Af)  Distribution 

Curve 


Fig.  3.4  Profiles  of  Stream  Function  Increment  Af 
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Fig.  3.5  Distribution  of  the  Derivatives  of  Af  (for  t=0.1)  and  Verification 
of  the  Satisfaction  of  the  Differential  EQuation  (...) 
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Fig.  3.8  Distribution  of  the  Shearing  Stress  at  the  Wall 


81 


Fig.  4.1  Boundary  Layer  Thickness  Over  a  Flat  Plate 

(Pr=l ,  y=l .4) 


Fig.  4.2  Boundary  Layer  Thickness  Over  a  Flat  Plate 

(Pr=l ,  y=l .4) 

(For  zero  heat  transfer  at  the  wall) 
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Transient  Contribution  to  the  Induced  Pressure  for  the  Case  of  No  Heat  Transfer  at  the  Wall 
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Fig.  4.3c  Transient  Contribution  to  the  Induced 
Pressure  for  the  Case  of  No  Heat  Transfer 
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Fig.  4.4a  Stepwise-decelerated  flow  Fig.  4.4b  Stepwise-accelerated  flow  (Ref.  [3]) 

Transient  Contribution  to  the  Induced  Pressure  for  the  Case  of  q  =  0.0 
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Fig.  4.4c  Transient  Contribution  to  the  Induced 

Pressure  for  the  Case  of  q  =  0.0 

w 


86 


»o 

O 

o 


o 


CM 

o 


CNOaO'O^J^OoO'O”^^ 

7  7  T  T  T  T  T 

e0lx(*d/dV) 


o 

m 


IX 


o 


Fig.  4.5a  Stepwise-decelerated  flow  Fig.  4.5b  Stepwise-accelerated  flow  (Ref.  [3]) 

Transient  Contribution  to  the  Induced  Pressure  for  the  Case  of  q  =  0.5 


Fig.  4.5c  Transient  Contribution  to  the  Induced 

Pressure  for  the  Case  of  q  =  0.5 
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Fig.  4.6a  Stepwise-decelerated  flow  Fig.  4.6b  Stepwise-accelerated  flow  (Ref.  [3]) 

Transient  Contribution  to  the  Induced  Pressure  for  the  Case  of  g  =  1.0 
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Fig.  4.6c  Transient  Contribution  to  the  Induced 

Pressure  for  the  Case  of  a  =  1.0 
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Fig.  4.7  Transient  Contribution  to  the  Induced 
Pressure  for  the  case  of  x.  =  3.5  and  M  =  10 
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Fig.  4.3  Transient  Contribution  to  the  Induced 
'Pressure  for  the  Case  of  x  =  3.5  and  =  20 


